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Stimulated Raman adiabatic passage (STIRAP) offers significant advantages for coherent population transfer 
between un- or weakly-coupled states and has the potential of realizing efficient quantum gate, qubit entangle¬ 
ment, and quantum information transfer. Here we report on the realization of STIRAP in a superconducting 
phase qutrit - a ladder-type system in which the ground state population is coherently transferred to the second- 
excited state via the dark state subspace. The result agrees well with the numerical simulation of the master 
equation, which further demonstrates that with the state-of-the-art superconducting qutrits the transfer efficiency 
readily exceeds 99% while keeping the population in the first-excited state below 1%. We show that population 
transfer via STIRAP is significantly more robust against variations of the experimental parameters compared to 
that via the conventional resonant n pulse method. Our work opens up a new venue for exploring STIRAP for 
quantum information processing using the superconducting artificial atoms. 

PACS numbers: 


Stimulated Raman adiabatic passage (STIRAP), which 
combines the processes of stimulated Raman scattering and 
adiabatic passage, is a powerful tool used for coherent popu¬ 
lation transfer (CPT) between un- or weakly-coupled quantum 
states^. It has been recognized as an important technique in 
quantum computing and circuit QED involving superconduct¬ 
ing qubits^ - — . For example, qubit rotations can be realized via 
STIRAP with the two computational states plus an auxiliary 
state forming a three-level A configuration^. A scheme for 
generating arbitrary rotation and entanglement in the three- 
level A-type flux qutrits is also proposed^, and the experimen¬ 
tal feasibility of realizing quantum information transfer and 
entanglement between qubits inside microwave cavities has 
been discussed^. Unlike the conventional resonant n pulse 
method STIRAP is known to be much more robust against 
variations in experimental parameters, such as the frequency, 
amplitude, and interaction time of microwave fields, and the 
environmental nois e 4 ^ 10 ? 11 . 

Recently, multi-level systems (qutrits or qudits) have 
found important applications in speeding up quantum gates^, 
realizing quantum algorithms^, simulating quantum sys¬ 
tems consisting of spins greater than l/2±£, implement¬ 
ing full quantum-state tomography^—, and testing quantum 
contextualityi^. Unlike the highly anharmonic A-type flux 
qutrits the phase and transmon qutrits have the ladder-type (E- 
type) three-level configuration which is weakly anharmonic. 
The dipole coupling between the ground state |0) and the 
second-excited state |2) in the phase qutrit is much weaker 
than those between the first-excited state |1) and the |0) state 
or the |2) state. In the case of the transmon qutrit the cou¬ 
pling is simply zero. This unique property makes it difficult 
to transfer population from |0) to |2) directly using a single n 
pulse tuned to their level spacing 6 l>2o- The usual solution is to 
use the high-power resonant two-photon process or to apply 
two successive n pulses transferring the population first from 


|0) to |1) and then from |1) to |2>I2J1. These methods often 
lead to a significant population in the middle level |1) result¬ 
ing in energy relaxation which degrades the transfer process. 
In contrast, STIRAP transfers the qutrit population directly 
from |0) to |2) via the dark state subspace without occupying 
the middle level |1). 

In this work, we report on the realization of STIRAP in a E- 
type superconducting phase qutrit—. As shown schematically 
in Fig. la the qutrit has a loop inductance L and a Josephson 
junction with critical current I c and capacitance C. The po¬ 
tential energy and quantized levels |0), 11), and |2) of the qutrit 
are illustrated in Fig. lb in which the frequencies to p ^ s of the 
pump and Stokes fields and their strength £l p , s (Rabi frequen¬ 
cies) are also indicated. Applying the rotating-wave approx¬ 
imation (RWA) in the double-rotating frame the Hamiltonian 
of the system can be written a s ^ L ^: 


H = 


o gp + g s e~ m 0 

g P + g s e i6t A p A(g p e iSl + gs) 

0 A{g p e~ lSt + g s ) A p + A s 


( 1 ) 


where the Planck constant h is set to unity, 6 = co p - to s , 
A p = coio-ojp and A s = (O 21 -to s are various detunings, g p , s are 
the qutrit-microwave couplings proportional to the amplitude 
of the pump and Stokes fields respectively, and A « V2 for the 
E-type configuration with weak anharmonicity. For 5 » Cl ps 
the fast-oscillating terms in equation (1) averages out to zero 
so the Hamiltonian becomes 



0 Q p /2 0 
top/2 A p L1J 2 
0 d s /2 A p + 


( 2 ) 


in which Cl p = 2 g p and = 2 Ag s . For the phase qutrit 
used here we have A ^ 1.45. Equation (0 is the well-known 
RWA Raman Hamiltonian^. In particular, when the system 
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FIG. 1: Superconducting phase qutrit and measurement pulse se¬ 
quences. (a) Schematic rf-SQUID type phase qutrit with Josephson 
critical current 7 C , shunt capacitance C, and loop inductance L. (b) 
Three bottom energy levels |0), |1), and |2) of the qutrit with related 
symbols indicated. Subscripts p and s refer to the pump and Stokes 
tones, respectively, (c) Counterintuitive pulse sequence with pre¬ 
ceding £l p for coherent population transfer from |0) to |2) without 
involving |1). (d) Conventional resonant n pulse sequence for suc¬ 
cessive |0) —» |1) —> |2) population transfers. 


satisfies the two-photon resonant condition: 


A p + A s — 0, 


(3) 


it has an eigenstate | D) = cos 0|O) - sin @|2), called the dark 
state, which corresponds to the eigenvalue of 6 = 0. Here 
tan 0(0 = £l p (t)/£l s (t). CPT from the ground state |0) to the 
second-excited state |2) without populating the first-excited 
state 11) can therefore be realized via STIRAP by initializing 
the qutrit in the ground state \0)^^ and then slowly increas¬ 
ing the ratio £l p (t)/Q. s (t) to infinity as long as the following 
condition s 1 i 2 i 24 i 25 


S » £lp, s , and 


rj 


Fl 2 p (t) + Q? s (t) dt > 107T (4) 


are satisfied so that the qutrit will stay in the dark state 
subspace spanned by ||0), |2)}. The first condition is required 
to reduce equation to equation (EJ) leading to the existence 
of the dark state solution while the second ensures adiabatic 
state following. 


Results 

Sample parameters and measurements. The sample used in 
this work is an aluminum phase qutrit, which is cooled down 
to T ~ 10 mK in an Oxford cryogen-free dilution refrigera¬ 
tor. The qutrit control and measurement circuit includes vari¬ 
ous filtering, attenuation, and amplification^. For the present 
experiment, we bias the rf-SQUID to have six energy levels 
in the upper potential well and use the lowest three levels as 
the qutrit states. The relevant transition frequencies are /io 
= oJio/2n - 5.555 GHz and /21 = (jl>2\I2ti = 5.393 GHz, 
and the relative anharmonicity is a = (/io - /2O//10 ~ 2.9%. 
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FIG. 2: Coherent population transfer via STIRAP in the super¬ 
conducting phase qutrit. (a) Stokes and pump microwave pulses 
Cl s (t) and Q p (t) shown in the overlapped region with the experimen¬ 
tal parameters Q 0 / 2n = 42.8+1.0 MHz and T d - 100 ns. (b) Level 
populations P 0 , P 1 , and P 2 versus time in the case of A p = A s = 0. A 
maximum experimental value of P 2 = 67%, limited by the low deco¬ 
herence times of the device, is reached. Experimental and calculated 
results are shown as symbols and lines, respectively. 


The measured energy relaxation times are 7/ 10 = I/T 10 = 353 
ns and T 21 = I/T21 = 196 ns, respectively, while the de¬ 
phasing time determined from Ramsey interference experi¬ 
ment is T*° = 124 ns. To realize STIRAP, a pair of bell¬ 
shaped counterintuitive microwave pulses with the Stokes 
pulse preceding the pump pulse, as illustrated in Fig. lc, are 
used. The pulses are defined by £l s (t) = £loF(t) cos[7r/(/)/2] 
and £l p (t) = D^F(t) sin[nf(t)/2] with F(t) = e ~^ /2Td ^ and 
f(t) = 1/(1 + e~ At/Td )^l. The pulse width (FWHM) is approx¬ 
imately 2T d , and its height is 0.967Qo, and £^(0) = Q p (0). 
Coherent population transfer. Figure 2a shows the two mi¬ 
crowave pulses defined by Qo/2tt = 42.8 ± 1.0 MHz and T d 
- 100 ns in their overlapping region. As t increases, £l s (t) and 
Q p (t) start to increase and decrease across t - 0 at which they 
are equal. The experimentally measured populations Pq, P 1 , 
and P2 versus time produced by this counterintuitive pulse se¬ 
quence in the resonant case A p = A s = 0 are plotted in Fig. 2b 
as symbols. We observe that as time evolves across t = 0 the 
population P2 (Pq) increases (decreases) rapidly, signifying 
the occurrence of STIRAP via the dark state of the supercon¬ 
ducting qutrit system. The experimentally achieved maximum 
P2, or the population transfer efficiency, is about 67% for the 
present sample under the resonant condition. Notice that in 
the entire region of t e [-100,100] ns, all of the characteristic 
features of the experimental data, in particular (i) P\ remain¬ 
ing significantly lower than P 2 for t > -50 ns, (ii) the decrease 
(increase) of P2 (Pq) after reaching the maximum (minimum) 
as well as the gradual rising of P 1 , are reproduced well by 
the numerical simulation (solid lines). The simulated tempo- 
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FIG. 3: Level population P 2 versus pump microwave detuning 
A p for four Stokes tone detunings. Bright resonances can be seen 
as the resonant condition equation (3) is satisfied (left-side peaks in 
each curves). A maximum value of P 2 = 72% is reached at the mi¬ 
crowave detunings A p = -A s = -20 MHz. The right-side peaks result 
from the two-photon process of the single pump microwave tone. 


ral profiles of the populations Pq, P i, and P 2 are obtained by 
solving the master equation p = ~(i/h)[H,p] + L(p ), using 
the measured qutrit parameters, where L(p) is the Liouvillean 
containing the relaxation and dephasing processes^. The nu¬ 
merical result also confirms that feature (ii) is due primarily to 
energy relaxation, while the maximum P 2 itself is mainly lim¬ 
ited by dephasing. Hence by increasing the qutrit decoherence 
times the undesirable decrease of P 2 and the rise of P 1 and Po 
in the relevant time scale can be suppressed and a higher P 2 
can be achieved (see below). 

In our experiment the conditions imposed by equation © 
are satisfied: 6/2n in the resonant case A p = A s = 0 is 
/10 ~ /21 = 162 MHz, which is approximately four times that 
of Do/2 tt, and it is easy to verify that the integrated pulse area 

JZ Tj£l 2 p (t) + £2*(t)dt* 327T is greater than IOjt. In addition 

to reducing decoherence the efficiency of the demonstrated 
STIRAP process can be improved by increasing the relatively 
small anharmonicity parameter a « 2.9% of the present sam¬ 
ple up to, say, 10% by optimizing device parameters of the E- 
type phase and transmon (or Xmon) qutrit s^^ . According to 
equation © greater anharmonicity allows the use of larger Do 
which would proportionally reduce the duration of the pump 
and Stokes pulses when the pulse area is kept unchanged to 
satisfy the adiabatic condition. Shorter pulses also reduce the 
negative effect of decoherence on the transfer efficiency. 

The STIRAP process is often identified in either the time 
domain or the frequency domain^. The latter is based on 
equation © which specifies the two-photon resonance con¬ 
dition. In Fig. 3, we show the level \2) population P 2 ver¬ 
sus the pump detuning A p for four Stokes tone detunings 
of As = 40,20,0, and -20 MHz, respectively. Bright reso¬ 
nance appears as the left-side peak in each curve when the 
two-photon resonant condition equation © is met. It is in¬ 
teresting to see that the maximum value of P 2 = 72% is 
reached at the microwave detunings of A p = -A s = -20 MHz, 


which is higher than the value achieved in the resonant case 
of A p = A s = 0 shown in Fig. 2b. In Fig. 3, the right-side 
peak in each curve is originated from the two-photon process 
excited by the single pump microwave tone. Compared to 
the left-side peaks, although the peak heights are comparable, 
they are much narrower, indicating that in practice it is less 
controllable using the two-photon process to perform coher¬ 
ent population transfer from state |0) to state \2). 

Efficiency and robustness. The above results demonstrate 
clearly CPT from the ground state |0) to the second-excited 
state \2) via STIRAP in the E-type superconducting qutrit. 
Compared to the usual high-power two-photon process or 
two non-overlapping successive resonant n pulse excitations 
shown in Fig. Id, which involve significant undesired popu¬ 
lation in the middle level |1) and require precise single pho¬ 
ton resonance and pulse areai2*±Z, CPT via STIRAP demon¬ 
strates just the opposite. First, in principle CPT between |0) 
and 12) can be accomplished without occupying the lossy mid¬ 
dle level |1). More importantly, the process is much more ro¬ 
bust against variations in the frequency, duration, and shape of 
the driving pulses^. In fact, in terms of the conditions equa¬ 
tion © and equation ©, we see from Fig. 3 that the single 
photon resonance condition is greatly relaxed. Although Q p s 
are limited by the system anharmonicity, their values, together 
with Td , still have much room for variations while maintaining 
the transfer efficiency. 

In Fig. 4a, we show the calculated results in the Do 
versus Td plane with contours indicating population transfer 
efficiency above 99% and 98%, respectively, using qutrits 
with decoherence times of T^ 10 = 35.3 yus, T^ 1 = 19.6 yus, and 
T™ = 12.4 //s, and the relative anharmonicity of a = 8%. 
The qutrit parameters in the ranges are now attainable with 
transmon 28 i 29 , Xmon—, flux—, and also possibly phase- 
type devices. In Fig. 4b, we show the level populations versus 
time (solid lines) for Do/27r = 100 MHz and Td = 50 ns, in 
which a nearly complete transfer above 99% from level |0) to 
level 12) is accomplished with the population in middle level 
|1) kept below 1%. These results indicate that the transfer 
efficiency of STIRAP is very insensitive to Do, which is 
limited by systems anharmonicity, and to Td, which should 
be much smaller than the decoherence time. The allowed 
variations of a few hundreds of MHz in Do and a few tens 
of ns in Td for keeping P 2 > 99%, for example, are in sharp 
contrast with the case of simple n pulse excitations. In fact, 
using a n pulse to transfer population from the ground state 
|0) to the first-excited state |1) alone, acceptable variations of 
the Rabi frequency D and pulse width T to keep the transfer 
efficiency above 99% can be estimated from the pulse area 
relation Dxf « (1 ± 0.06)7r. Thus if we use D/27T = 50 
MHz, the pulse width variation must be less than ± 0.6 ns. 
More strict condition is required when successive excitation 
from the first-excited state |1) to the second-excited state |2) 
is considered. From these we see that the extreme robustness 
of the STIRAP process is very advantageous and should be 
useful in various applications such as realizing efficient qubit 
rotation, entanglement, and quantum information transfer in 
various superconducting qubit and qutrit systems. 
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FIG. 4: Calculated results demonstrating the robustness of the 
STIRAP process, (a) Results in the G 0 versus T d plane with con¬ 
tours for population transfer efficiency above 99% and 98% indi¬ 
cated. The same sample parameters as those in Fig. 2b are used in 
the calculation but the decoherence times are increased to T\° = 35.3 
ps, 7\ 21 = 19.6 ps, and T*° = 12.4 ps and the relative anharmonicity 
to a = 8%. (b) The temporal profiles of microwave pulses (dashed 
lines) and level populations (solid lines) in the case of Clo/2n =100 
MHz and T d = 50 ns. 


Discussion 

We have experimentally demonstrated coherent population 
transfer between two weakly-coupled states, |0) and |2), of a 
superconducting phase qutrit having the E-type ladder con¬ 
figuration via STIRAP. The qutrit had a small relative anhar¬ 
monicity of 2.9%, and moderate decoherence times of T*° 
= 124 ns, Tj 10 = 353 ns, and T^ 1 = 196 ns, respectively. We 
demonstrated that by applying a pair of counterintuitive mi¬ 
crowave pulses in which the Stokes tone precedes the pump 
tone, coherent population transfer from |0) to \2) with a 72% 
efficiency can be achieved with a much smaller population in 
the first-excited state |1). Using the measured qutrit parame¬ 
ters, including decoherence times, we simulated the STIRAP 
process by numerically solving the master equation. The re¬ 
sult agrees well with the experimental data. We showed that 
by increasing the decoherence times of the qutrits to the order 
of a few tens of microseconds, currently attainable experimen¬ 
tally, the transfer efficiency can be increased to greater than 
99% while keeping the population of the first-excited state be¬ 
low 1%. 

We have also shown that coherent population transfer 
via STIRAP is much more robust against variations of the 


experimental parameters, including the amplitude, detun¬ 
ing, and time duration of the microwave fields, and the 
environmental noise over the conventional methods such as 
using high-power two-photon excitation and two resonant n 
pulses tuned to oqo and oj 2 \, respectively. Therefore STIRAP 
is advantageous for achieving robust coherent population 
transfer in the ladder-type superconducting artificial atoms 
that play increasingly important roles in various fields 
ranging from fundamental physics to quantum information 
processing. Furthermore, the method can be readily extended 
to the A-type systems such as the superconducting flux 
qutrits, in which the initial and target states locate in different 
potential wells representing circulating currents in opposite 
directions. Our work paves the way for further progress in 
these directions. 


Methods 

Determination of level populations. The qutrit level popula¬ 
tions P[{t) (i = 0,1,2) at a given time t are determined using 
two carefully calibrated nanosecond-scale measurement flux 
pulses A and B, which reduce the potential barrier to two dif¬ 
ferent levels so that tunneling probabilities p A and p B in each 
case are measured. Pulse A leads to low tunneling probabil¬ 
ity Pq(~ 5%) for state |0), high tunneling probability p A for 
state |1) and, of course, even higher tunneling probability p A 
for state \2 ). Pulse B results in a slightly deeper potential well 
than pulse A does so that p B ^ 0, p B « 5%, and a much larger 
p B for state |0), |1), and |2) respectively. Denoting the density 
operator of the qutrit as p, we have 

P A ’ B = PoPo’ B + Pipf B + PiP A 2 B , (5) 

where Pf = pa , and p A,B can be found from the experimen¬ 
tally determined tunneling probabilities pi of the ith energy 
level for given amplitudes of pulses A and BA&. Combining 
the normalization condition 


Trp = P 0 + Pi + P 2 = 1 , (6) 


we obtain 


Pi = l(Pj ~ Pk)P A + (P A k ~ P A )P B + PjPk ~ P A kPj 1 ID , (7) 

where i - 0, 1,2 with {/, j , k} in circulative order like {0,1,2}, 
{1,2,0}, and {2,0,1}, and D is the determinant 


D = 


i »a 


1 1 

/i js 

P i 

Pi 


(8) 


Hence the level populations Pq, Pi, and P 2 can be obtained by 
measuring p A and p B . 

Numerical simulations. We numerically calculate the level 
populations P 0 (t) = poo(0, P\if) = pn(t), and P 2 (t) = p 22 (t) at 
any given time by solving the master equation 


p = ~[H,p\ + L(p), (9) 

where p is the system’s 3x3 density matrix, H is the Hamil¬ 
tonian given by equation (1), and L(p ) is the Liouvillean con¬ 
taining various relaxation and dephasing processes. Since ex¬ 
perimentally the pump and Stokes microwaves are not cor¬ 
related, we introduce a phase difference 0 between the two 
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microwaves in the actual calculations^. In this case, the wave approximation leads to a Hamiltonian in the following 
double-rotating reference frame is described by the operator form: 

U = |0)(0| + |1><1|^ + |2><2 and the rotating- 

I 


H = 


0 gp+gse-W-® 0 

g P +gsettgpeW-*) + g s ] 

0 Mgpe-**-** + g s ] A p +A,. 


( 10 ) 


where the Liouvillean operator in equation (9) is given by^: 


J 


L(p) 


loyi 

(Tio + / 




v-*- iu ■ /i 0 )pio 2ri 0 pii - 2r 2 ip 2 2 (rio + r 2 i + )pi2 

(r 2 i + r 2 0 )P 1 2 3 4 5 ° (Ho + r 2 i + 72i)P2i 2 r 2 ip 22 


(ii) 


We use the parameters Tiq = 2.83xl0 6 7 8 9 sec -1 , T 2 \ = 5.10xl0 6 
sec -1 , and = 8.06xl0 6 sec -1 measured directly from ex¬ 
periment, while from the measured dc Stark shift of the two- 
photon spectral lines we estimate y * 0 ~ 2y^ 0 and y^i ~ 7\o- 
In our calculations p(p 0) is obtained by solving equation (9) 
using the fourth-order Runge-Kutta method. Since the phase 
difference 0 of the two microwaves in our experiment is ran¬ 


dom, we average the result over 0 and finally arrive at: 

1 r 2n 

pit) = 2tt J 0 P(*>0)#- (12) 
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